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ABSTRACT 

Modelica is a new object-oriented multi-domain modeling lan- 
guage based on algebraic and differential equations. In this paper 
we present an environment that integrates different phases of the 
Modelica development lifecycle. This is achieved by using the 
Mathematica environment and its structured documents, "note- 
books". Simulation models are represented in the form of struc- 
tured documents, which integrate source code, documentation and 
code transformation specifications, as well as providing control 
over simulation and result visualization. 

Import and export of Modelica code between internal; struc- 
tured and external textual representation is supported. Mathematica 
is an interpreted language, which is suitable as a scripting language 
for controlling simulation and visualization. Mathematica also sup- 
ports symbolic transformations on equations and algebraic expres- 
sions, which is useful in building mathematical models. » 

BACKGROUND 

Integrated simulation environments are advantageous in or- 
der to work effectively and flexibly with simulations. Users 
prepare and run simulations as well as investigate simulation 
results. Several auxiliary activities accompany simulation 
experiments: requirements are specified, models are de- 
signed, documentation is associated with appropriate places 
in the models, input and output data as well as possible con- 
straints on such data are documented and stored together 
with the simulation model. The User should be able to repro- 
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duce experimental results. Therefore input data and parts of 
output data as well as the experimenter's notes should be 
stored for future analysis. 

Traditionally, simulation and accompanying activities 
have been expressed using heterogeneous media and tools: 

• a simulation model is traditionally designed on paper 
using traditional mathematical notation; 

• simulation programs are written in a low-level 
progranmiing language and stored on text files; 

• input and output data (if stored at all) are saved in 
proprietary formats needed for particular applications 
and numerical libraries; 

• documentation is written on paper or in separate files 
that are riot integrated with the program files; 

• the graphical results are printed on paper or saved using 
proprietary formats. 

When the result of the research and experiments, such as a 
scientific paper, is written, the usernoiimally gathers togeth- 
er input data, algorithms, output data and its visualizations as 
well as notes and descriptions. One of the major problem in 
siinulation development environments is that gathering cor- 
rect versions of all these components from various files and 
formats is difficult and error-prone. 

USING MATHEMATICA NOTEBOOKS 

Our approach to the integration problem is based on the 
Mathematica (Wolfram 1997) environment arid its program- 
mable notebooks. Every notebook corresponds to one docu- 
ment (one file) and contains a tree structure of cells. A cell 
can include other cells and/or arbitrary text or graphics. In 
particular a cell can include a code fragment or a graph with 
computatioiial results. The hierarchy of cells corresponds to 
the traditional hierarchy of chapters, sections, and para- 
graphs used in textual documents and available in advanced 
word processors like FrameMaker or MSWord. The text can 
be written in different styles, including certain font families, ^ 
bold, italic, color and size as well as mathematical typeset- 
ting. 

The contents of cells can be 

• parts of models (a formal description that can be used for 
verification, compilation and execution of some 
simulation model); 

• text/documentation (used as comments to executable, 
formal model specifications); 

• dialogue forms for specification and modification of 
input data; 

• result tables (the results can be immediately represented 
in table form); 
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• graphical result representation (with 2D vector and 
raster graphics as well as 3D vector and surface 
graphics); 

• 2D graphs that are used for various model structure 
visualizations: 

• . class diagrams 

• variable dependency diagrams 

• data structure diagrams 

THE MODELIC A LANGUAGE 

The language called Modelica (Modelica 1998) for hierar- 
chical physical modeling has been developed in an interna- 
tional effort. It is an object-oriented language (Elmqvist, 
1997, Fritzson and Engelson, 1998) for modeling of physical 
systems. The language unifies and generalizes previous ob- 
ject-oriented modeling languages. Modelica is intended to 
become a de facto standard. It offers three important fea- 
tures: 1) non-causal modeling based on differential and alge- 
braic equations; 2) multidomain modeling capability, i.e. it 
is possible to combine electrical, mechanical, thermodynam- 
ic, hydraulic, etc. model components within the same appli- 
cation model; 3) a general type system that unifies object- 
orientation, multiple inheritance, and templates within a sin- 
gle class construct. 

Modelica- niodels are built from classes: Like in other 
object-oriented languages, a class contains variables, i.e. 
class attributes representinjg data, The main difference 
compared to traditional object-oriented languages is that 
instead of functions (methods) the programmer uses 
equations to specify behavior. Equations can be written 
explicitly,, like a=b. or can be inherited from other classes. 
Equations can also be specified by the connect statement. 

Modelica model of lunar landing 

As an introduction to Modelica we will present a model of a 
rocket landing on the moon surface. 

There are two components (Rocket and Planet) used in 
the Landing class. The class Rocket e;mbodies floating 
point (Real) variables, simulation-time constants 
(parameter) and equations for vertical motion of the rocket 
above the planet surface. The motion is influenced by 
gravitational force and rocket motor force (Thrust). Rocket 
mass varies since the rocket loses fuel mass proportional to 
the amount of thrust from, the rocket motor. Velocity and 
acceleration are time-derivatives of height and velocity as 
specified by the der (...) construct. 

The class x^anding includes apollo (an instance of 
class Rocket) and moori (an instance of Planet, with 
specified planet mass and radius). . . 

Dynamic systems are described by models where 



behavior evolves as a function of time, i.e. all the variables 
in the model evolve as functions of time. Modelica has a 
predefined variable time which steps forward during system 
simulation. The simulation models the rocket behaviour 
from time =0 until the rocket touches the ground; 

Equations describe the thrust function, a step function 
with an initial thrust force level fl during the time interval 
[0. . thrustDecreaseTime] and a second thrust level f2 
during the time, interval [thrustDecreaseTime.. 
thrustEndTime]. The Step functions can conveniendy be 
expressed using Modelica conditional expressions (the if - 
then-eXse construct). 

The gravitational field is expressed as: 



(h(t) + radiuspi3„gt) 

The Lunar landing model expressed in Modelica . 

class Rocket "generic rocket class" 
Real height (start=59000) ; 
Real velocity(start=-2000); 
Real mass (start=1000) ; 

Real acceleration; . ^ 
Real Thrust; , 
Real gravity; 

parameter Real massLossRate=0 . 000277 ; 
e<3uatlon 

Thrust-nvas^*gravity=mass* acceleration; 
dor (mass) =-massLossRate*abs (Thrust) ; 
der (height ) =velocity; 
der (velocity) =acceleration; 
end Rocket; 

class Planet "generic planet class" 

parameter Real mass; 

parameter Real radius; 
end Planet; 

class Landing "landing of a rocket onto a 
planet" 

parameter Real f orcel=36000 ; 
parameter Real force2=1500; 
parameter Real thrustEndTime=210; 
parameter Real thrustDecreaseTime=43 ; 
parameter Real G=6.672e-ll; 
Rocket apollo ; . 
Planet moon (mass=7 . 3 82e22 , 
. radius=1.738e6) ; 

e(7uatlon 

apollo . Thrust= 
if ( tijne< thrustDecreaseTime) then forcel 
else if (time< thrustEndTime) then force2 

- else 0 ; 

apollo.gravity=(G*moon.mass) / 

(apollo . height +m9on . radius) ^^^2 ; 
end Landing; 
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Note that Model ica equations do not specify which variables 
are inputs and which are outputs. Thus, the causality of equa- 
tions-based models is unspecified. This is fixed only when 
the equation system is solved. 

Simulation Semantics 

Classes, instances, and equations are translated by the Mod- 
el ica compiler into a flat set of differential-algebraic equa- 
tions, constants and variables. The initial values for tiine=0, 
specified by the ( s tart= . . . ) construct can be taken from 
' the model definition. A simulation engine finds a numerical 
solution of the system of ordinary differential equations. All 
model variables are functions of (modeled) time during the 
simulation. 

The simulation engine is a fairly complex piece of 
software that can be implemented in many different ways. In 
each case this engine sorts equations (Elmqvist, 1997), finds 
algebraic loops, checks consistency of ordinary differential 
equations (ODE), controls an ODE solver, and computes 
requested variable values from the solution found. 

Existing Environment and Need for MathModelica 

There is an environment for the Modelica language and en- 
gine for Modelica simulations, based on the Dymola system 
(Dymola 1998). The code, documentation, input data, graph- 
ical and numerical results are represented in different, heter- 
ogeneous formats, which is a disadvantage of the system. 
Furthermore, the user is not able to extend the system and in- 
troduce new software components. Additionally, the user 
cannot perform user-specified symbolic (algebraic) rhanipu- 
lations with the formulae used in the code. 

For these reasons, we are developing the extensible 
integrated MathModelica environment which is partly based 
on the Mathematica system and notebooks. 

A specific feature of Mathematica is that models are 
normally not written as free formatted text. Instead, 
Mathematica expressions (terms) are used. These can be 
conveniently written in a tree-like prefix form, or entered 
using standard mathematical notation. Every term is a 
number; an identifier or a form such as: 

head [termj, . . - # term^l 

In order to satisfy this requirement, we designed the new 
MathModelica language^. Note that MathModelica has the 
same abstract syntax and the same semantics as Modelica, 
but a different concrete syntax. This means that essentially 
the same language constructs are written differently, as illus- 
trated below. 

The MathModelica language uses some Mathematica 
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notation, such as: 

termj/ . - - ; term^, {terznj, ... - , teirm^} , 

termi term2, termi^^term2, term', termi"^ term2^ 

and arbitrary arithmetic expressions composed from terms. 

The ' (tick) means time-derivative. The ^ (backtick) 
character in MathModelica corresponds to . (dot) notation in 
Modelica. - 

Lunar Landing Example in MathModelica 

Class [Rocket ; 
Declare [ 

Real [start->59000,unit->-in'' ] heights- 
Real [start->-2 000,unit->'*iri/s'' ] velocity; 
Real[start-> 1000, unit-> "kg" ] mass; 
Real acceleration; 
Real Thrust; 
Real gravity; 

Parameter Real inassLossRate=0 . 000277 ; 

]; 

EQuatlon [ 

Thrust-inass*gravity==inass*acceleration; 
mags ' ==-massLQSsRate*Ab^ [Thrust] ; 
height ' ==velocity; 
velocity ' ==acceleration; 

] 

J'- 

Class [Planet; -**********************- 
Declare [ ' 

Parameter Real mass; 
Parameter Real radius; 

3 

BeginClass [Landing] ; 
Declare [ 

Parameter Real forcei=36000 ; 

Parameter Real force2=1500; . 

Parameter Real thrustEndTime=210 ; 

Parameter Real thrustpecreaseTim.e=43 ; - 

Parameter Real G=6 . 672*10^-11 ; 

Rocket apollo; • 

Planet [mass->7 .382*10^22, 

radius->1.738*10''6] moon; 

]; 



Eqiiation [ 

apollo * Thrust == 
Whicli [ tijne<thrustDecreaseTime , f orcel , 
tijne<thrustEndTime, . f orce2 , 

True , . 0 ] ; 
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apollo *gravity= = 

■ G • moon 'mass ' 

(apol lo ' height + moon * radius)^ 
EndClass [Landing] ; 

Expressions used in MathModelica are "dynamic". This 
means that they can be created as result of symbolic and al- 
gebraic transformation performed by Mathematica. These 
transformations can be integration, differentiation, expan- 
sions of series, simplification, etc. 

Formal Syntax of MathModelica 

The formal concrete syntax of MathModelica is described by 
a BNF gramhiar^ as shown below, where keywords are all 
marked with bold, other terminals are courier (if they are 
written as here) or italic (generic tokens, like integers or 
strings): 

program: := class* 

class: := shortClass | longClass 

shortClass : : = Class [ciassjiame; 

(ext I decls (equations )* ]; 

\yhen the short class definition is evaluated in Mathematica, 
the class description is added to the model. 

longClass : : =BeginClass Iclassname} ; 

(ext I decls | equations | comment )* 
EndClass iclassname'} s 

The constructs BeginClass [] and EndClass [] are bal- 
anced. Class definition can be split, and written in several 
cells. The EndClass {name} ; adds (when evaluated), defi- 
nition of class name to the syntax tree. These constructs al- 
low the user to insert documentation into class code 
(Figure 1). 

A class can extend (inherit) other classes as in Java: 

ext ::= Extends iclassname*} s 
decls : := Declare [ onedecl* 1; 

equations ::= EquationC oneequation* ]? 
onedecl : := [Parameter] type 
[ [ attr *] ] [ [idx *3 ] 
varname [ = value ];[ "short comment" ; ] 

Note that "pariameters" are simulation time constants and 
their value can be stated in the simulation model.. 

idx ::= integer (used for array declaration) 
attr : : = . 
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start-> value; (bound value of variable for time=0) 
|.unit->" string" (unit used for user interface only) 
. I varname- >value ; (specifies value for component) 

comment : : = text written in a separate notebook cell 

type : := Real | Integer | classname 

oneequation : : = expr expr j 

expr : : = varname^ | expr + expr 

I expr- expr | expr/ expr |' expr* expr 
I JbuiJt injunction (expr) | expr [[expr]] 
[ expr ^varname | number j expr®^^ 



Fig. 1. Notebook fragment with MathModelica class declaration 
and comments. 

In Figure 1 the symbols in the right margin denote different 
kinds of cells and the hierarchy constructed from the cells. 
The bracket □ means a cell with code (an executable com- 
mand); the bracket with the horizontal bar 3 means a cell 
with text (a fragment of documentation). A larger □ -brack- 
et marks a cell that contains eight smaller cells. A hierarchy 
of cells of unlimited depth can be easily created. 

THE EXPERIMENTAL ENVIRONMENT 

MathModelica is both a language and an environment. 

The Mathematica built-in functions and MathModelica 
environment functions can be used in order to manipulate the 
model. Naturally, the most important operation using the 
model is simulation. This can be performed in different 
w^ays: 

• translation of the model to Mathematica, and simulation 
using the Mathematica built-in ODE solver; 

• translation to Modelica. and using the Dymola-based 
environment for simulation; 

• translation to C-f+ using the MathCode system 
(MathCode 1998, Fritzson 1997), a Mathematica to G+-i- 
compiler. 
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ln[23]:=BeginClass [Lemdlng] ; ^ 




This class describes the landing process and evaluates | 
the height of the rocket and other dynamic parameters _| 




.ln[24]:=Declare [Parametex- Real £o2rcelB3 6000; ] 




This force variable is needed... ^ 




ln[25]:=Declare [Parazaeter Real £orce2ol500; ] ; ~j 




This force variable is needed too... 3 




In[26]:=EndClass [r.aiiding] ; 





TRANSLATION FROM MATHMODELICA TO 
MATHEMATICA 

An order to convert a MathModelica model to Mathematica 
the function ModelicaToMathematica [ ] is called. Vari- 
ables and functions declared in the model then become inter- 
actively accessible by the user working in the Mathematica 
environment. 

Classes, variables and equations in the MathModelica 
model are converted to a flat list of variables and equations. 
First, the constants are set: 

setparam[ ] = (apollo ^massLossRate = 0.000277; 
thrustDecreaseTiine=43 ; 
thrust EndTiine=2 10; 
moon *mass=7 . 3 82*10 ""22; 
moon * radius =1 . 738*10^6; 
forcel=3 60d0; 
force2=1500; 
G=6. 672*10^-11) ; 

The Mathematica differential equation solver (NDSolve) re- 
quires that only differential equations are specified in the 
equation list, so algebraic equations have to be. handled sep- 
arately. 

There is an equation (here named eql) with the variable 
apollo* acceleration the translator has to extract: 

eql: = 

apollo * Thrust [t] - apollo *mass [ t] * 
. apollo ' gravity [t] =^ 

apoilo'mass [t] *apollo * acceleration [t] ; 

To extract the variable apollo * acceleration the alge- 
braic equation solver is called: 

Solve[eql,apollo * acceleration [t] ] ; 

The Solve function expresses the variable apollo* ac- 
celeration [t] in terms of other variables. A new func- 
tion definition for the acceleration is created: 

apollo * acceleration] [t_] : = - ( ( 

apollo * gravity [t] * apollo *mass [ t ] - 
apollo *Thrust [t] ) / apollo *mass [ t ] ) ; 

Two other functions are already given in the MathModelica 
code: 

apollo * gravity [t_] : = (G*inoon * mass) / 

(apollo * height [t] +moon * radius) ^2; 

aE>ollo*Thrust[t_] : = 

If[t < thrustDecreaseTime ,forcel, 
IfCt < thrustEndTime, force2, 0]]; 



eqs={ \ ■ 
apollo *mass' [t]==- apollo *massLossRate* 

TU^sC apollo * Thrust [t] ] , 
apollo * height ' [ t] ==apollo ^velocity [t] , 
apollo ^velocity ' [t] ==apollo *acceleration[ t] 

}; .... 

Initial conditions are necessary for the solution: 

initcond={ 

apollo ^height [01 ==59000, 
apollo Velocity [0] ==-20a0 , 
apollo'mass tO]==1000}; 

Solutions are required for these functions: 

vars={apollo * height [t] , apollo 'velocity [t] , 
apollo * mass [t 3 ) 

When the function Model icasimula te [timejna^] is 
called, initial conditions, equations, and functions are passed 
to an ordinary differential equation solver which is a Mathe- 
matica built-in function: 

res=NDSolve [ Join [eqs , initcond] , var , 
{t, 0, timejnax)^ 

As a result. Mathematica creates polynomial approximations 
for all functions we want to find, known as interpolate-func- 
tions. 

Result visualization 
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Fig. 2, The. apollo * height value for time in [O^timejuaJ 
(timeinax=500) 

The functions computed by NDSolve can be visualized 
graphically. For instance, the height of the rocket can be vi- 
sualized (see Figure 2) by the call 

ModelicaPlot [ apollo * height [t] , 

"Lxinar landing: Altitude"]; 



Now, three ordinary differential equations are specified: 



other Scripting Utilities 

Three scripting utilities - ModelicaToMathemaitica [ ] , 
ModelicaSiinulate[] and ModelicaPlot [ ] have been 
illustrated above. 

iinportModelica[] accepts textual fraditional 
Modelica code, parses it and converts to MathModelica 
code. ExportModelica[3 perfornis the opposite 
translation. It can be used with a Modelica implementation, 
available in the Dymola environment. 

Input Data Specification 

Another script conunand, ModelicainData [ 1 allows se- 
lecting parameter values and specifying boundary conditions 
interactively. The dialogue window is created automatically 
from specification of data structureis (Engelson and Fritzson, 
1996) in the Modelica code (see Figure 3). Variable values 
can be edited and new simulation runs can produce results 
immediately. In the future this can be used for computational 
steering of Modelica simulations. 



Fig^ 3. A Mathematica dialogue box for setting inpqt variable 
values of a MathModelica simulation. The boxes apollo and 
moon represent cell instances with several components. 

CONCLUSION 

The MathModelica is an extension of the Modelica lan- 
guage targeted for work within the Mathematica environ- 
ment. This language is object-oriented (the programs consist 
of collections of classes). The language is equation based: 
instead of traditional functions and procedures we use non- 
causal equatioiis which specify algebraic and differential re- 
lations between numerical variables. Input-output causality 
is not specified, and therefore these equations can be used in 



multiple ways. 

The environment integrates most activities needed in 
simulation design and use: documentation, modeling 
(coding), symbolic processing and transformation of 
formulas, input and output data visualization. This advanced 
programming environment can be applied in various 
simulation applications. 
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